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We theoretically evaluate the squeezed joint operators produced in a single optical parametric os- 
cillator which generates quadripartite entangled outputs, as demonstrated experimentally by Pysher 
et al. [T][Phys. Rev. Lett. 107, 030505 (2011)]. Using a linearized fluctuation analysis we calculate 
the squeezing of the joint quadrature operators below threshold for a range of local oscillator phases 
and frequencies. These results add to the existing theoretical understanding of this potentially 
important system. 
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I. INTRODUCTION 



II. SYSTEM AND METHOD 



Cluster states are highly entangled multipartite states 
that have attracted much interest since they were pro- 
posed by Raussendorf and Briegel [2j as a resource state 
with potential applications for one-way quantum com- 
puting. However, as with conventional approaches, issues 
of scalability arise [3j [4] . Recent experimental advances 
have been made by Pysher et al. [1 in the generation 
of a continuous- variable (CV) quadripartite cluster state 
from a single optical parametric oscillator (OPO), based 
on a proposal by Menicucci et al. [5 j. In this approach 
quantum registers are encoded in the quadratures of the 
optical frequency comb produced by the spectrum of the 
OPO. The experimental scheme was shown to be scal- 
able and capable of simultaneously generating 15 quadri- 
partite entangled cluster states. In previous work with 
others [6j 0, we verified the presence of quadripartite 
entanglement in such a system by demonstrating viola- 
tions of the van Loock-Furusawa (VLF) criteria [8 j. We 
also calculated the squeezed joint quadrature operators 
in the temporal domain to determine whether the state 
produced could be classified as a cluster state. 



In this Brief Report, we investigate the squeezed joint 
quadrature operators in the frequency domain as this is 
the form more amenable to comparison with experiment. 
A spectral analysis of the system is required as the tem- 
poral domain solutions give the squeezing integrated over 
all frequencies and are not useful when we wish to know 
for which frequencies the system will be most useful. By 
verifying that the squeezed joint quadrature operators 
approach very low noise levels at particular frequencies, 
we determine that the scheme approaches an ideal cluster 
state. This Brief Report is structured as follows. Section 
II provides details of the theoretical model and method. 
Section III sets out the results and discussion. A sum- 
mary of our findings is given in Sec. IV. 



The system considered here is a pumped optical cav- 
ity containing a x^ non-linear crystal, as investigated 
theoretically in Ref. [7]. We extend our analysis to in- 
corporate the results of a full spectral analysis in the 
frequency domain. The optical cavity is pumped by two 
field modes which give rise to four low-frequency output 
modes at frequencies CJ3, CJ4, CJ5, ujq. Mode 1 is pumped 
at a particular frequency and polarization such that it 
produces modes 3 and 6, as well as modes 4 and 5. Mode 
2 is pumped such that it gives rise to modes 5 and 6. 

The master equation for this system is [9 



dp 

dt 



D 

Him + H pump , pj + ^ 'Yi'Di [p] (1) 



where the Hint is the interaction Hamiltonian given by 



Hint = ih[xiaia\a\ + xi^i^l + ^^aja^] + h.c, (2) 

with the Xi representing the effective nonlinearities and 
cti and a\ denoting the bosonic annihilation and creation 
operators for the intra-cavity modes at frequencies c^. 
We note here that we are using the full interaction Hamil- 
tonian rather than making a classical undepleted pump 
approximation leading to linear Heisenberg equations as 
used by Pysher et al. [1 (see also the Supplemental Ma- 
terial), which has the advantage that our formalism may 
also be used above the oscillation threshhold |6j [7] . In 
a rotating frame, the pumping Hamiltonian, H pU mp is 
given by 



H n 



z 



(3) 



where €i are the classical laser amplitudes, the ji repre- 
sent the cavity losses at each frequency and T>i[p] is the 
Lindblad superoperator [9] describing zero-temperature 
Markovian damping. 

In order to investigate the squeezed joint quadrature 
operators, and in turn determine whether the entangled 
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output beams emerging from the cavity approximate a 
four mode cluster state, we begin by defining the quadra- 
ture field operators for each mode and arbitrary phase 
angle, as 



and for the low frequency fields we find 



cue 



Yi = Xi(0 + 7r/2). (4) — = Xi«i«5"-74«4 



The squeezed joint quadrature operators are then the 
eigenvectors of the system found by solving the Heisen- 
berg equations of motion from the interaction Hamilto- 
nian in the undepleted pump approximation. These are, 

01 = (-c 1 X 3 + c 1 X 4 -X 5 +X 6 )e- C2r , 

2 = (-c 2 X 3 -c 2 X4 + X 5 +X 6 )e- Cir , 

3 = (c 1 Y 3 + c 1 n + Y 6 + Y 6 )e- ( » r , 

04 = (c 2 %-c 2 Y A -%^%)e-^ r , (5) 

where c t = (y/E - l)/2, c 2 = (y/b + l)/2 and r is the 
squeezing parameter. The common eigenstate of these 
operators tends towards a quadripartite entangled cluster 
state when r — >• 00. The graph for the corresponding 
cluster state is depicted in Fig. [I] 
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FIG. 1: (color online) Graph of a weighted square-cluster 
state corresponding to the squeezed operators in Eq. |5| . The 
vertices represent the intracavity modes, the edges represent 
the coupling between these and fractions are the weightings. 



A. The positive-P representation equations 

We use the positive-P representation [10] equations for 
the system described in Sec. [TTJ These are derived using 
the standard approach [12] whereby the master equa- 
tion is mapped onto a Fokker-Planck equation for the 
positive-P function using a correspondence between two 
independent stochastic fields ol{ and ol\ and the mode op- 
erators di and aj, respectively. The resulting equations 
can then be written as a set of c-number stochastic differ- 
ential equations. This method allows us to calculate clas- 
sical stochastic averages which are equivalent to quantum 
mechanical normally-ordered operator moments. 

The positive-P equations for the high frequency fields 
are 



da 1 

dOL2 



e 1 - Xl«4«5 - Xi«3«6 - 7i^i, 



C2 - X2«5«6 - 72^2, 



(6) 
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Xia^a^ -73^3 



Xi«i 
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(779(f) +%o(t)), 



M*)+"7e(*)), 



Xl^iaJ + X2«2«6" _ 75^5 



X2«2 



Xl^ia^ + X2«2«5" _ 76^6 



Xi^i 



(jlb(t) - irj 6 (i)) 



(779 (t) - ivio(t)) 



X2&2 



(m(t)-iV2(t)). 



Since the phase-space is doubled there are also the equa- 
tions found by swapping with a^~ and rji(t) with 
r]i+2(t)- The rji(t) are real, independent, Gaussian noise 
terms with correlations rji(t) = and rji(t)r}j(t f ) = 
5 ZJ 5(t-t'). 



B. Linearized Fluctuation Analysis 

We now use Eqs. ([6| and ^ as the starting point 
for a linearized fluctuation analysis [9 . Linearizing 
these around the classical steady-state solutions, we find 
the evolution equations for the fluctuations. Specifi- 
cally, we neglect the noise terms in Eq. ^ and set 
di = OLi + 5di, where ai are the steady-state val- 
ues and 5ai represent the fluctuations. To first or- 
der the equations of motion for the fluctuations, Sex = 
[Sai^Saf, 5a 2: Sa^, . . • , SaQ, Sa^] 7 ', can be written as [7 



dSa = -AScxdt + BdW , 



(8) 



where B is the noise matrix of Eq. ^ with the classical 
steady-state values inserted, dW is a vector of Wiener in- 
crements [12] and A is the drift matrix with the classical 
steady-state values inserted. 

Provided that the eigenvalues of the drift matrix A 
have no negative real part, the system is well described 
as a multi- variable Ornstein-Uhlenbeck process [TT] . for 
which the intracavity spectral correlation matrix is found 
as a matrix product which is then related to the measur- 
able output fluctuation spectra by the standard input- 
output relations [T3] , 



III. RESULTS AND DISCUSSION 

We now calculate the spectral variances of the squeezed 
joint operators and confirm that these approach zero over 
a range of phases and frequencies. We consider the sym- 
metric system where the two pumping inputs, ei 5 2, are 



(7) 
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equal, the two nonlinear it ies, Xi,2 are equal, and the cav- 
ity losses, 7i,..,6 = 7, are equal. As shown in Ref. 0, the 
system has an oscillation threshold and for the parame- 
ters used here the critical pump amplitude is e c = 61.8. 
The parameters used in the following calculations are 
7 = 1, x — 0-01 and e = 0.94e c , giving a pump rate ap- 
proximately 12% below threshold in intensity. For these 
parameters, this was found to yield the best squeezing. 

In Fig. [2] we plot the minimum variance of the squeezed 
joint operators, V(Oj), as a function of the phase at a 
fixed frequency. We find the frequency at which the min- 
imum variance is obtained and plot all four variances at 
this corresponding frequency. Specifically, Fig. [2] repre- 
sents the maximum squeezing achieved by the operators 
V(0\$) but not ^(02,4). The frequency that yields the 
maximum squeezing corresponds to uj = 0.357, corre- 
sponding to 30 MHz with a cavity loss rate 80 MHz. We 
find that V{0\) = V(Os)^ as indicated by the (green) 
solid line, and V{02) = V{0^)^ as indicated by the (blue) 
dashed line. The corresponding flat lines indicate the 
coherent state level for the variances of the particular 
quadrature combinations, the dashed line being at 7.24 
and the solid line at 2.76. This result shows squeezing 
well below the coherent state level and identifies a region 
where the squeezed joint operators approach zero for a 
range of phases centered around = and multiples of 
7r. This provides a theoretical quantum treatment of the 
experimental results obtained in Ref. pQ and provides 
further evidence that an ideal continuous variable four 
mode cluster state can be closely approximated by their 
single OPO scheme. 
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FIG. 2: (color online) Plot of the minimum variances of the 
squeezed joint operators as a function of phase at a fixed 
frequency of uj = O.357. 

Figure [3] depicts the same quantity shown in Fig. [2] as 
a function of phase, but on a logarithmic plot. This plot 
compares well to the experimental results presented in 
Ref. pQ with the width of the squeezing region becom- 
ing smaller as the degree of squeezing increases. The 
squeezed joint operators are scaled relative to the rele- 



vant coherent state level of zero. 
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FIG. 3: (color online) Logarithmic plot of the minimum vari- 
ances of the squeezed joint operators as a function of phase 
at a fixed frequency of uj = 0.357- 

We also investigate the behavior of the squeezed joint 
operators over a range of frequencies. Figures [4] and [5] 
show plots of the minimum variances as a function of 
both frequency and phase. The range of parameters over 
which the entangled state can be considered as approx- 
imately a four mode cluster state can be seen in the 
smooth region surrounding = and where V(Oj) — >• 0. 
In both Fig. [4] and Fig. [5] the functions are cropped at 
V(Oj) = 8 to improve visualization. 




FIG. 4: (color online) Plot of the minimum variances of the 
squeezed joint operators, V(Oi) and V(Os), as a function of 
phase and frequency. 

Finally, we consider the regime of optimal squeezing 
for the joint operators. Figure [6] shows the minimum of 
the sum of the variances of both squeezed joint operator 
pairs as a function of phase on a logarithmic scale. Al- 
though this compromises the overall degree of squeezing 
obtained, this regime may be of use in experiments as the 
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IV. CONCLUSIONS 




FIG. 5: (color online) Plot of the minimum variances of the 
squeezed joint operators, V(02) and V(Oa)^ as a function of 
phase and frequency. 

corresponding frequency of uj = 0.237 yields the best pos- 
sible degree of squeezing for the combined quadratures. 



In conclusion, we have demonstrated that the re- 
quirements of a CV four mode cluster state are closely 
satisfied in the single OPO scheme proposed in Ref. [5] 
and experimentally realized in Ref. [1 . This is verified 
by confirming that the squeezed joint quadrature opera- 
tors for the system approach zero for a range of phases 
and frequencies. These findings are in good agreement 
with the experimental results of Pysher et al and 
extend current theoretical results pQ into the spectral 
domain. This work builds on our previous findings in 
Ref. [7] wherein CV quadripartite entanglement was 
demonstrated in the system and the squeezed joint 
quadrature operators were calculated in the temporal 
domain. Using our approach, it is a simple matter to 
extend the analysis to parameter regimes above the 
oscillation threshhold, which may have advantages when 
photon flux is important. 
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